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Abstract — Stochastic advection diffusion equation (SADE) with multiplicative stochastic input is a practical mathematical 
model for different physical phenomena. In this paper, SADE will be studied using two spectral stochastic techniques. The 
first is the Wiener chaos expansion (WCE) technique and the second is the Wiener-Hermite expansion with perturbation 
(WHEP) technique. These techniques convert the SADE into a system of deterministic partial differential equation (DPDE) 
that can be solved using a deterministic numerical method which is suitable for the periodic boundary conditions. 
Convergence analysis is discussed and some of the second order moments are compared. The numerical results demonstrate 
the efficiency of both techniques. The WCE technique is more accuracy than the WHEP technique. The diffusion and 
advection coefficient and the intensity of Gaussian white noise play important roles in the SADE solution. The study shows 
that the WCE technique is more practical to get the closed form mean solution while the WHEP technique gets the mean 
solution in the form of an infinite series. 

Keywords — Advection diffusion equation , Multiplicative random input , Wiener-Hermite expansion, Wiener-Chaos 
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I. Introduction 

Numerical models have received a great attention in sciences and engineering in the recent years for modeling the differential 
equations. These models work on reducing the cost and time of computation especially for the physical phenomena that 
contain uncertain input. This phenomena can be studied by converting it into mathematical models of stochastic differential 
equations (SDE) and we can use the numerical methods [1,2] for overcoming these problems. 

Numerical methods have been developed for simulating SDE such as moment equations, probability density method [3], etc. 
These methods are complicated in solving the nonlinear SDEs so the spectral decomposition techniques have received much 
attention in the recent years. The spectral decomposition technique was first suggested by the great mathematician Norbert 
Wiener [4]. Wiener constructed an orthonormal random basis for expanding homogeneous chaos depending on white noise, 
and used it to study problems in statistical mechanics [5]. The Hermite polynomial has been used to obtain the solution of 
SDE. Mecham et al. [6] suggested the Wiener-Hermite expansion to study turbulence solution of Burger equation. In 
nonlinear stochastic differential equations, there exist always difficulties of solving the resultant set of deterministic integro - 
differential equations. The deterministic integro-differential equations got from the applications of a set of comprehensive 
averages on the stochastic integro-differential equation obtained after the direct application of WHE. Many authors 
introduced different methods to face these obstacles. Among them, the WHEP technique [7] was introduced using the 
perturbation technique to solve perturbed nonlinear problems. M. El-Tawil and his co-workers [7-11] used the WHE together 
with the perturbation theory (WHEP technique) to solve a perturbed nonlinear stochastic differential equation. The WHEP 
technique is generalized to handle n th order polynomial nonlinearities, general order of WHE and general number of 
corrections [8]. 

Cameron and martin [12] was developed a more explicit and intuitive formulation for the Hermite polynomial, which was 
called the Wiener-Chaos expansion. Their development is based on an explicit discretization of the white noise process 
through its Fourier expansion. This approach is much easier to understand and more convenient to use, and hence replaced 
Wiener’s original formulation. Fourier chaos expansion has become a useful tool in stochastic analysis involving Brownian 
motion [13]. Rozovskii et al. [14-16] derived Wiener chaos propagator equations for several important the stochastic partial 
differential equations (SPDEs) driven by Brownian motion forcing. Lototsky et al. [17, 18] proposed a new numerical 
method for solving the Zakai equation based on its Wiener chaos expansion. Using Fourier-Hermite expansion for modeling 
non-Gaussian processes is also investigated [19, 20]. Babuska et al. [21], Schwab et al. [22] and Keese el al. [23] developed 
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and generalized Ghanem's approach for solving stochastic elliptic equations. Xiu and Karniadakis [24] generalized the 
Hermite polynomial expansion and used it to study flow-structure interactions. Zhang et al. [25] combined moment 
perturbation method with polynomial chaos expansion, and used it to study the saturation flows in heterogeneous porous 
media. 

The main goal of this paper is to use two stochastic spectral techniques, WCE and WHEP for solving the stochastic advection 
diffusion equation with multiplicative white noise and periodic boundary conditions. The two techniques convert SPDE into 
a system of DPDE. The DPDE can be solved using a proposed Eigen function expansion in both cases of WCE and WHEP 
techniques. The results will be studied through the mean and variance solutions. 

This paper is organized as follows: 

The formulation of the SADE is outlined in section 2. The WCE technique is explained in section 3. In section 4, the 
algorithm of WHEP technique is introduced. In section 5 and 6, we apply the WCE and WHEP techniques respectively; the 
convergence analysis of the WCE is studied. The proposed method for solving the resulting DPDEs; the numerical solutions 
of the WCE and WHEP techniques are introduced in section 7. The comparison and discussion of the results of the two 
techniques are in section 8. Finally, the conclusions are given in Section 9. 

II. Stochastic Advection Diffusion Equation 

SADE represents the transporting that occurs in fluids through the combination of advection process and diffusion process 
.Consider the SADE in Stratonovich sense with multiplicative stochastic force described by white noise and periodic 
boundary conditions as [27]: 

^2 ^ 

du(x, t)= a —- u(x , t)dt + g — u(x , t) o dW ., t > 0 , 

dx 2 dx (2.1) 

u( 0, x) - sin(x), x g (0, 2n), 
which can be written in the Ito sense as: 

d 2 d 

du(x,t ) = y — -u(x,t)dt + <7 — u(x,t)dW. t> 0 , 

dx 2 dx (2.2) 

u( 0, x) = sin(x), x g (0,2;r), 

where u(x,t) represents the concentration of mass transfer; JLI> 0 represents the diffusion coefficient; CJ > 0 

2 

CJ 

y = ju + and yv t is the one dimensional Brownian motion. 

This model has exact solution [27] as: 

u(x , t) = e~ ^ sin(x + cjW( t)). (2.3) 

The first order moment E[u(x , t )] and the second order moment E[u(x , t )] 2 are: 


E[u(x,t)\ -e r sin(v) 


E[u(x 9 t)] 2 = e 


-2 /dt 


l l 

- e 

V2 2 


-2 at 


cos(2x) 


(2.4) 


So the solution admits the Wiener chaos expansion because of finite second moment. 

III. Wiener Chaos Expansion Technique 

For SPDE with random force in form of Brownian motion (BM) and for fixed time T > 0 , we consider an orthnormal basis 
in Hilbert space L 2 ([0,7"]) to be the trigonometric functions [16]: 
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m 1 ( t ) = J , m i (t) = J— cos 


(i - X)nt 


i > 1, 0 < t < T. 


V T ) 

Also define the standard Gaussian random variables C t (GRVs) and the Brownian motion W t = W (t) as follows: 


1 00 1 

£ =\m i {s)dw s , = 

0 0 

The convergence in mean square sense [26] in the interval [0, T ] is: 

E[W{t)-^\ mi {T)dTf <^, 0 <t<T 


(3.1) 


=1 0 


According to [12], the solution function u(x, t) can be expanded as follows: 

u{x,t) = Y J u a (,x,t)r a 00, u a = E[uT a (£)] = E[uT a (t)], 


(3.2) 


where the multi-indices 3 = {a = (a i ,i > 1), a i e {0,l,2,...};|a| = T a , < 00 } and T a (C) — ek (Ci) where H n (x ) 

i =1 t =1 

is the normalized n th order Hermite polynomial and T a are called Wick polynomials of order OC . Also 7^(7)is the Wick 
polynomials filtered by the a-algebra F f W , Recall that T a (t) is a martingale and satisfies the differential equation [26]: 


dT a (0 = Z dW > ’ 


(3.3) 


where multi-index OC i (j) = 


(Xj, J*l, 

a t _i, j=i. 


Wick polynomials form a complete orthonormal basis in the Hilbet space. The expectation of two wich polynomials T a , 7^ 
is E[TTJ = S , 

L a (3 J a ,p 


v (K + N)\ 

Truncating the expansion (3.2) up to polynomial of order N and using only K GRVs retains ^ 7 ^—jj coefficients [26]. 


The truncated multi-indices will be N = {a = (cq,.. .OC K ), a i e {0,1,2,... }',\oc\ = < N}. Then the truncated 

i =1 

WCE can be denoted as: 

u KN (xj)= Yj u a( xd K(C)- (3.4) 

ae3 K,N 

The mean and variance of the truncate solution function U K N (x, t ) are computed as: 

E[u KN (x,t;co)] = u 0 (x,t\ 

2 

var [u K N (x,t;a>)]= Z K(-*V)| • 
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IV. Wiener Hermite Expansion 

As a consequence of the completeness of the Wiener-Hermite set [28] any arbitrary stochastic process can be expanded in 
terms of the Wiener-Hermite polynomial set H w (t 19 ...,t n ) . This expansion converges to the original stochastic process with 
probability one. 

The stochastic solution process u(x, t\ w) can be expanded as, [7]: 


u(x , t ; co) = u (0) (x, t) + ^ J u (k) (x, t\t x 

k =1 R k 


t k )H {k \t v ...,t k )dT k , 


(4.1) 


where dr k = dt l dt 2 ....dt k and J is a ^-dimensional integral over the disposable variables t l ,t 2 ,...,t k ■ The functional 

R k 

H (n) (t v t 2 ,...t n ) is the n th order Wiener-Hermite time independent functional and i > 0 are the 

deterministic kernel of the WHE. 


dW 


The Wiener-Hermite functionals // (77 form a complete set with H (0) =land H a \t) =- is the white noise. The H 


r(n) 


dt 


functions are statistically orthnormal, i.e. 

£[/f <0 ] = 0 V/ > 1, 

e\ 


~H (n H U) ~\ = 0 Vi9ty. 


(4.2) 


The solution will be practically truncated with (m +1) terms and the expectation and variance of the truncated solution will 
be: 


E[u(x , t\(D)\ = u (0) (a, t ), 

m 2 

var[M(x,f;ry)] = ^/:![ [u (k) (x,t,t v ...,t k )} dr k . 

k=\ R i 

In the nonlinear SPDE or the multiplicative SPDE, it is difficult to solve the deterministic differential-integral equation 
system of the kernels results from the application of the WHE. This difficulty is due to the resulting system is a coupled 
differential-integral system and we can overcome this by using the perturbation technique. 

In the perturbation technique the solution is a power series of small parameter (7 . Set of simple equations are expanded as 
[29]: 

NC 

u (k) = ^ & u\ k) ,k> 0, (4.3) 

*•=0 

where NC is the number of corrections. For ffl order WHE, the statistical properties of the relatively solution will be 
calculated as: 


NC 

E[u(x,t;a>)\ = yV'uf 0 (x, t ), 

i =0 

m f NC 

var[u(x,f;ft>)] = V&! [ Vcr'u^Cx,?;^ 

k 4 V i =0 


\ 2 

t k ) dz k . 


) 


V. Application of WCE 

To get the WCE of SADE (2.2); consider the differential form: 


(4.4) 
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d[uT ] = du.T + u.dT + du.dT , 

L a J a a a 7 

using equations (2.2) and (3.3) to get: 


d[uT] ■■ 


^ fi 2 

r 5? 


d ( /— 

u(x , t)dt + a — u(x, t)dW t T a +u 2^'\J a i m i (0^- ( .) 


+ 


^ 8 2 

r s* 


o c v—i l 

u(x, t)dt + cr — u(x , f)dW f . 2^ *\j a i m i (0 T- (i) dW t 


(5.1) 


Taking the expectation for both sides of (5.1), the terms involving Ito integrals will disappear since they are mean zero. 
Then we will get, 


d[u(x,t)] = 


\ 


y^TU a (x,t)dt 

OX 


+ CT 


U 


V i> 1 


(5.2) 


. . . . _ f sin(v), a- 0, 

with the initial condition U a (X, UJ — < q 


Theorem: 

The SADE (2.1) which has exact solution u(x,t) and the truncated solution u K N (x,t) ; the estimated error will take the 
form: 

-jut s&T (<T T) N+1 


E\u(xj)-u KN (xJ)\<Ce ^ + 


K (N + l)\ 


Proof: 


according to [27], since u(x , t ) = 6{x + oW ( t ), 0? and VE(0 ==- 7 = Ci + 

/=1 0 1-2 

w( 0 =w*+w s . 


Vir (; - \)Kt 

- sin(-), let 


(/ — I )7T 


W = — 1 = C +y, 4 ' ^ sin(f-—-—) and IV = Y ^ sin(f-—-—) where W K and W s are orthogonal. After 

Vr « '(i-l)a- r ' (i—l)*r T 

that we can expand the solution #(x + oW K + <jW s ,t) by Taylor's series with respect to W K and W s respectively, to get: 


u(x, t )= o(x, t)+ y, 


(aW K ) m d" 


: 0(x, t)-\ - 0(x + <jW k + rh , t)cjw s + 




K aiitx r r $ j 

V'+l ^N+l 


m=1 m! fix'" fix 

Truncating the solution with respect to order N and K GRVs to get: 


(N + l)! dx 


N +1 


0(x + r/ 2 ,t). 


u 


K,N 


, N . ^ i . ow K ) m dr Q . , 
(x, 0=< 9 (x, 0+y —f— —#(x, 0, 

m=t ™ fix 


2 (T 

then we have E u(x , t)-u„ AT (x, t) < C, oE W„ + - E W. 

K ’ N 1 Y ^ +1 / AT , 1M 7 


(iV + 1)! 


where C = sup 

ft x X 


fix 


0(x,t) 
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2 T 


1 


Ai-l)jct s 2 T 


then E\W;\<— Y--sin 2 (--—) < 

1 s| n 2 T rfK 

E | W 2N+2 1 = (E | W 2 |) N+1 (2N +1)!! < (E | W 2 |) A ' +1 (2N +1)!!, 

where ( 2N +1)!! = ( 2N + 1)(2 N -1).1, 

then, 

gT (a T) n+1 


E u(x,t)-u KN (x,t) <Ce ^ 


- + - 




K (7V + 1)! (53) 

This means that, the error between the approximate solution and the exact solution decays by increasing the number K of 
GRVs, the order of polynomial chaos N and diffusion coefficient /U . Also decreasing the advection coefficient G and the 

time interval T increases the convergence between the exact and the WCE approximation. 

VI. Application of WHE 

The first two terms in expansion (4.1) are the Gaussian part of the solution. This part is not sufficient for the accurate 
solution of SADE. The second order WHE is applied to the SADE (2.2) to get: 


— u ( 0 \x,t)+ f— u {l) (x,t\t x )H {X) (t x )dt x + f f— u ( 2 \xJ;t x J 2 )H ( 2 \t x J 2 )dt x dt 
p+ j p+ J J Pf 


ed 


d 


dt 


r 


dt 


dt 


' d 2 


‘ ' d 2 


dx 
crn(t) 


—^u {0 \x,t) + j—^u {1) (x,t;t l )H {1) (t l )dt 1 +|| —^■u {2 \x,f,t l ,t 2 )EI {2) (t 1 ,t 2 )dt l dt, 


dx 


+ 


— u {0) (x,t)+ [ — u m (x,t;t l )H m (t l )dt l + [ [ — u' 2> (x,t;t t ,t 2 )H' 2> (t l ,t 2 )dt l dt 2 
dx {dx J J Av 


0 

(°) 


dx 


0 0 

u {V) (x, 0) = sin(x), u a) (x, 0; t x ) = u {2) (x, 0; t x , t 2 ) = 0. 


( 6 . 1 ) 


Multiplying both sides of (6.1) with H i0 \H a) (t x ) and H^ Z) (t x ,t 2 ) respectively and taking the expectation to get: 


( 2 )/ 


— u (0) (x, t) = y — -u (0) (x, t) + G — u (l) (x, t\ t). 
dt dx dx 


— u (1) (x, t\t x ) = y —- u (l) (x, tUi) + G — u (0) (x, t)S(t -t x ) + 2 g — u {2) (x, t\ t,t x ). 
dt dx dx dx 


( 6 . 2 ) 


(6.3) 


d d 2 \ d \ d 

— u {2) (x, t’,t,,t 0 ) = y —- u {2) (x, t\t , l ) -f— cr — u {l) (x, t;t)S(t -t 0 )-\— g — u {l) (x, t; L )8(t -1 ,). 

/"s . L Z ' ^ 2 1Z 1 Z Z I 


dt 


dx 


2 dx 


2 dx 


(6.4) 


The deterministic systems appear when applying WCE and WHE are coupled integro-differential system of equations that 
are not easy to solve. In the following section, we will suggest a numerical technique to solve the deterministic systems (5.2) 
and ( 6 .2-6.4) with periodic boundary conditions. 

VII. The numerical technique using the principal of Eigen function expansion 

We introduce a numerical technique for solving the DPDE resulting from the application of the WCE and WHEP techniques. 
Here we implement the usual Eigen function [30] to be suitable for solving the resulting system of differential equations with 
periodic boundary conditions. The solution u(x , t ) can be written as an Eigen function expansion as: 
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<(*,*) = 2 


n =0 


t 

C 1 n e~ an 2 t +je- an 2 (t - s) F n 1 (s)ds 


where, 


Cg=jjf(x)dx, 


C\ = — J/(x)cos(nx)dx, 


Cl = — J f(x)sin(nx)dx, 


cos(nx ) + ^ 


ft=l 


t 

C 2 <T"" 2 ' + j e^ a " 2 (l ^ > F i l (s) ds 


sin(nx), 


1 L 

Fj = —J F(x,t)dx, 


Fl = — | F(x, t ) cos(nx) dx, 


F n 2 = — | F(x, t ) sin(nx) dx. 


(7.1) 


(7.2) 


The proof of formulae (7.1) is explained in Appendix A. 

The system of propagators (5.2) can be re-written using the Eigen function expansion (7.1) as: 


t (*,o= jr 


n =0 


i 

C 1 n e~ an2 ' +fe~ an 2 (, ~ s) F n ‘(s)ds 


cos(nx) + ^ 


n =1 


i 

Cl e an2t + Je - 0 " 2 (> ~ S) F 2 (s ) * 


For 1^1 = Owe get: 

0 0 2 

— u(x , 0 = 7 —7 w(x, 0? (X 0) = sin(x), 

dt dx 2 

which results in: 

u a (x, t ) = e~ n sin(x). 

For |Of| = 1, <x x = 1, a i =0we get: 

d d 2 d 

— U(X,t) = r -r-T M ( X ’ 0 + OT,1 1 (0 M U|=0 ( X ’ ? )’ U a ( X > 0 ) = 

dt ox ox 1 

which results in: 

(x, 0 - e~ yt cos(x). 

\JT 

For |6)r| = 1, a t = 1, a x = 0, i > 1 we get: 


sin(ftx). 


(7.3) 


(7.4) 


(7.5) 


0 0 2 0 

— u(x,t) = y —- w(x, t) + crm. (0 — u, , 0 (x, 0, O, 0) = 0, 
dt dx dx 

which results in: 

/2 

u a (x, t ) = — sin(£;zt) cos(x), k = (i — 1) / 7\ 




kslf 


(7.6) 
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For oc = 2, a { = 2, a i =0, i > 1 we get: 


— u(x,t) = r —^u(x,t ) + crmj(?) — w, , (x, t), 

dt ox ox 11 ’ 1 


u a (x, 0) = 0, 


which results in: 


u a (x,t) 


~<J t 


4it 


e rr sin(x). 


(7.7) 


For \a\ = 2, a l = 0, a t = 2, i > 1 we get: 

9 9 2 (— d 

— u(x,t) = y —- u(x , 0 + V2 am (0 — w, , =1 a =1 (x, 0, (x, 0) = 0, 

dt dx dx 

which results in: 


u a (x,t) 


-V 2 o - 2 _ 

Tn 2 k 2 * 


sin 2 (£;z7) sin(x). 


(7.8) 


For a = 2, a } = 1, cr =1 , / >1 we get: 


d d 2 d d 

— u(x,t) = y —- w(x, 0 + am l (t) — w, , x a =1 (x, 0 + crm. (t) — m, , t a =1 (x, f), (x, 0) = 0, 

dt dx dx dx 


which results in: 


u a (x,t) 


-y/2 cr 2 t , 

- e 

T n k 


sin (knt) sin(x). 


(7.9) 


For \a\ — 2, a i —l, a [ — 1, i, l > 1 we get: 

8 8 2 8 .8 

— u(x,t) = r — u(x,t) + crm i (t) — u w , la (x,t) + cnnj(t)—u M=ia=l (x,t), u a (x,0 )=0; i,j> 1, ; < ./', 
of ox Ox ox 

which results in: 


u a (x,t) 


- 2 -<? sin(/c7rr)sin( tt.s' ) sin(x), k 2 = j-l/T. 

T 7t k k 2 


(7.10) 


Then we have: 


u 0 =e r ‘ sin(x). 


(7 _ 


m h=i "Tf 


e yt cos(x) 


sfl . ' 

t H-sin(£;z7) 

7lk 


, fc = (/-i)/r, 


W M =2 “ ' 


£ 7 sin(x) 


t 2 727 ^ . 2 . 72 _ : _ 2 


—H-sin(A;;zT) + 

V2 n & n kk. 


sinC&^zt) sin(& 2 ;zT) + 


n 2 k 2 


sin (k nt) 


,k 2 = j-l/T. 

The WHE differential equations (6.2), (6.3) and (6.4) can be solved using perturbation technique [7]. We can use the 
perturbation technique about advection coefficient a combining with the Eigen function expansion (7.1) to calculate the 
first, second, third and fourth corrections. 
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Compare the coefficients of cr° : 


8 


8^ 

u°(x,t)- J u— T u°(x,t) = 0, Uq(x,0) = si 
dt ox 

o 


sin(x). 


to get ill (X t) = e sin(x), 

d d 2 

-u l 0 (x,t;t l )- j u—ru l 0 (x,t;t l ) = 0, u l 0 (x,0\t l ) = 0, 
dt ox 

to get u {) (x,t\t ] ) = 0 , 

d 2 d 2 2 2 

—u 0 (x,t;t l ,t 2 )- J u— T u 0 (x,t;t l ,t 2 ) = 0, u 0 (x,0;t v t 2 ) = 0, 
ot ox 

to get w,j (x,t;t l ,t 2 ) = 0. 

Compare the coefficients of <7* : 

d 2 d 

ul (x, t)- ju —j Mj° (x,t) = — u\ (x, t ; f), Mj° (x, 0) = 0, 
dx dx 


d 

- LI , 

Sr 


to getMj°(x,f) =0, 
d d 2 d d 

— ul(x,t;t l )- J u—ru{(x,t,t l ) = S(t-t l )—ul(x,t) + 2—ul(x,t,t,t l ), m^x.O;^) = 0, 
ot ox ox ox 

to get u\ (x, t\t x ) = e~ M ‘ cos(x), 

d 2 Id Id 

u 2 (x,t\t x ,t 2 ) - jU—ruf (x,f,t v t 2 ) = -5{t -t 2 )—ulixj;^) + - S(t -t x )—ul(x,t;t 2 ), 
dx 2 dx 2 dx 

2 tv rw : t = n 


dt 

uf (x, 0; t x , t 2 ) = 0, 
to get uf (x,t;t r t 2 ) = 0. 

2 

Compare the coefficients of <J : 


d 0 d 2 o 1 d 2 o d i 

— u 2 (x,t)-ju—ru 2 (x,t) = -—ru 0 (x,t) + — u l (x,t;t), 
dt dx 2 dx dx 


u° 2 (x, 0) = 0, 


_3 

to get u 2 (x,t) - — t e~^ sin(x). 


d d 2 \ d 2 d d 

— u i 2 (x,t;t l )- J u — u i 2 (x,t;t l ) = - — u 1 0 (x,t;t l ) + S(t-t l ) — ul(x,t) + 2 — uf(x,t,t,t l ] 
dt dx 2 dx dx dx 


u\ (x, 0; ) = 0, 
to get u 2 (x,f,t\) =0, 


(7.11) 


(7.12) 


(7.13) 


(7.14) 


(7.15) 


(7.16) 


(7.17) 


(7.18) 
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d 8 2 Id 2 18 

u 2 (x,f,t x ,t 2 )- ju—^ul(x,t',t l ,t 2 ) = ——^ul(x,t,t v t 2 ) + — 8(t-t 2 )—u\(x,f,t l ) 


dt ' dx 2 


to get u 2 (x,t;t l ,t 2 ) = -e M sin(x). 


Compare the coefficients of <7 : 


2dx 2 


dx 


1 d 

+—S{t—t x )— u\{x,t\t 2 ), ul(x,0;t x ,t 2 ) = 0, 

2 dx 


dt 


d d 2 Id 2 d 

u° 3 (x, t)- ju u 2 ( x, t) = - j Mj° (x, t ) + — u\ (x, t ; t), u° 3 (x, 0 ) = 0 , 


dx 1 


2 dx 


dx 


to get ll\ (x, t) = 0, 


d i. . d 2 j. , 1 d 1 j. ... . d o. . . d 2 . 

u 3 (x,t',t x )- ju—^u 3 (x,t\t j) = — —^u l (x,f;fj) + o(f-f 1 ) —m 2 (x,/) + 2 — u 2 (x,t;t,t x ), i 


dt ‘ dx 2 


2 5x 


5x 


5x 


5 3 

to get u 3 (x,t;t l ) = ——te M ‘ cos(x) - — t x e ^ cos(x), 

d d 2 \ d 2 i d 

—ul(x,t;t l ,t 2 )-ju—rH 3 (x,t;t l ,t 2 ) = -— i: uf(x,t;t l ,t 2 ) + -S(t-t 2 )—u l 2 (x,t;t l ), 
dt dx 2 dx 2 dx 

1 d 

+ — 8(t — t x ) — u\(x,t\t 2 ), u 3 (x,0;t l ,t 2 ) = 0, 

2 dx 

to get u 2 (x,t',t l ,t 2 ) = 0. 

Compare the coefficients of <7 4 : 

d d 2 Id 2 d 

— u° 4 (x, t)~JU —T u° 4 (x, t) = - —r u 2 (x, t) + — u\ (x, t\ t ), U° 4 (x, 0 ) = 0 , 
dt dx 2 dx dx 


to get u° 4 (x, t) - — 1 2 e sin(x), 
8 


d i d 2 i 1 5 2 ! 5 0 / x ^ 2 / \ 

—M 4 (xa;f 1 )-//^-yM 4 (x,f;f 1 ) = -^- u 2 (x,t;t l ) + o(t-t x ) —u 3 (x,t) + 2—u 3 (x,t,t,t x ). 


dt ' dx 2 


2 dx 


dx 


dx 


to get i 4 (x, 0 ;f,) = 0 , u\(x,t;t l ) = 0. 


d 2 , . d 2 2 , , 1 C 2 , , L, . a !, 

— U 4 {x,t,t v t 2 )- JU — U 4 (x,f,t l ,t 2 ) = -—u 2 {X,t;t 1 ,t 2 ) + -d(t-t 2 )—u 3 (x,t;t 1 ) 


dt 11 '5x 2 


2 ax 2 


2 ^ dx 

. 2 , 


1 ^ 

+ -£(7-^) — u 3 (x,t;t 2 ), u 4 (x,0;t l ,t 2 )=0, 

2 dx 


to get iC 4 ( x > * ’ h ’ *2 ) = ~ sin(jc) + 2 t 2 e jL " sin(x) + 2 t x e Ml sin(x). 


Then we have 


(7.19) 

(7.20) 

''(x, 0 ;/|) = 0 , 

(7.21) 

(7.22) 

(7.23) 

(7.24) 

(7.25) 
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x 19 

u (0) (x,t) = e~ M> sin(x)(l — a 2 t-\ - cr 4 t 2 ), 

2 8 

( ^ 3 'N 

w (1) (*,0 = <re _/ " cos(jc) 1 — <r 2 (r + —f t ) , 

V 2 5 j 

u {2) (x,t ) = <j 2 e~ M ‘ sin(x)[-l + 0.5 cr 2 (t + 4fj + 4r 2 )). 

Then the mean and the variance are: 

3 19 

2?[w(jc,f)] = u {0) (x,t) = e~ Mt sin(x)(l - <j 2 t-\ - cr 4 ^ 2 ), 

2 8 

t t t 

var (w(x,0) = | [w (1) (x,f;f 1 )] 2 ^ 1 + 2 J J [M (2) (id;^d 2 )] 2 ^ 25 

0 0 0 


2-2 Lit 

— g e 


COS 2 (jt) f f - — CT V + — <rV 

l 2 4 J 


+ 2ct 4 £ 2// ' sin z (x) 




5 23 83 44 

— err 4- — crV 

2 24 


(7.26) 


(7.27) 


The WHEP technique gives the mean and the variance of the solution in the form of an infinite series in G . Since G 
represents the perturbation parameter in the WHEP technique, decreasing this parameter gives a good convergence. 

VIII. Numerical Results 

In order to examine the efficiency of the proposed methods, comparisons between the approximate solutions and the exact 
solution of the SADE are simulated through the following figures. We also introduce some discussion about the effect of the 
diffusion coefficient and advection coefficient. 

We take T = 5 and 50 Gaussian random variable. As shown in Figs (1-22), we can note, the approximate solution of the two 
methods and the exact solution are in satisfactory agreement with each other under some convergence conditions. The 
convergence of the WCE is enhanced by increasing the number of Gaussian random variables K . The convergence of the 
WHEP is enhanced by decreasing the perturbation parameter G which represents the advection and also white noise 
coefficient. 


T=5;cf=0.01 

Exact mean 



T=5;a=0.01 

Exactvariance 




FIG 1. The exact mean at x = n I 3 for different FIG 2. The exact variance at x = n I 3 for different 


values of jU # 
T=5;a=0.01 

Exact mean 



values of jU # 


T=5;<7=0.01 


Exactvariance 
0.00003 



1 2 3 4 5 6 

Fig 4. The exact variance at t = 0.5 for different 


of//. 


values of // _ 
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Exact mean 


T=5;/i=2 


FIG 5. The exact mean at x - n 13 for different 
values of* 7 . 


T=5;p=2 


Exactvariance 



1 



FIG 6. The exact variance at x = 7r 1 3 for different 
values of cr 


First, we examine the effect of the diffusion coefficient jU and the advection coefficient G on the exact solution with respect 
to the convergence rate (5.3) in Figs (1 - 6). Studying those figures, we can note that, for fixed point X, the decaying 

exponential function e ,ut affects the solution. The mean solution and the variance solution decay faster with time at large 
values of diffusion coefficient jil . The larger the diffusion, the vanished faster the variance over time is. For fixed point t , 
the effect of the sinusoidal function appears on the behavior of the solution. As the diffusion coefficient jU increases, the 
variance decreases with time. 


We investigate the effect of the advection coefficient G on the solution of the mass transfer in Figs (5, 6). Examining those 
figures we found that, as the coefficient G increases, the solution decreases and vanishes with the time. As the coefficient G 
decreases, the variance also decreases. So choosing large value of fi and small value of G is appropriate for good 
convergence of the approximate solution. 


T=5;GRV=50;cF=0.01;pr=2 

mean 



T=5;GRV=50;a=0.01;//=2 

variance 



_ Exadvarianoe 

— IstWCEvarianoe 

- - IstVVHEvariance 


Fig 7. Second order means response for the exact, 
WCE and WHEP. Comparison between the three 
means at x = 7r / 3 . 


Fig 8. First order response variances for the exact, 
WCE and WHEP. Comparison between the three 
variances at x = n / 3. 


T=5;GRV=50;cf=Q.01;p=2 

mean 



Exactamean 

WCEmean 

WHEmean 


T=5;GRV=50;cf=0.01,/j=2 

variance 



— Exactv aria nee 

15tV'/C E va ria n c 

1 atVVHEvarianci 


FIG 9. Second order means response for the exact, 
WCE and WHEP. Comparison between the three 
means at t = 0.5 . 


FIG 10. First order response variances for the exact, 
WCE and WHEP. Comparison between the three 
variances at t = 0.5 . 
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T=5;GRV=50;a=0.01;//=2 

variance 





— Exa ctvaria n ce 

— 2nd WC Evaria n ce 

— - 2ndWH Evaria noe 


Fig 11. Second order response variances for the 
exact, WCE and WHEP. Comparison between the 
three variances at x = n / 3 . 


T=5;GRV=50;<7=0.01;ps2 

variance 

7.* 

6 . 

5.* 

4.x 
3.x 
2 . 

1.x 





Exactvariance 

2 ndWCEvaria nee 

2ndWH Evaria n oe 


Fig 12. Second order response variances for the 
exact, WCE and WHEP. Comparison between the 
three variances at t = 0.5 . 


The comparison between the exact solution and both approximate solutions of the WCE and WHEP techniques are shown in 
Figs (7-11). Examining those figures elucidate a satisfactory agreement between the exact solution and the approximate 
solutions either in mean or variance. 

The effect of increasing the magnitude of the advection term GU x is examined in Figs (13-16). We can note that increasing 

G makes a deviation between the exact solution and the approximate solutions especially in the case of using the WHEP 
technique. This deviation is a result of using the WHEP technique which gives a mean solution and variance in the form of 
an infinite series in G and t . 


T=5;GRV=5-0;cf=1 ;jj=2 

mean 



Exactmean 
WCEmea n 
WHEmean 


T=5;GRV=50;«cf= 1 ;pr=2 
variance 



— 

— Exa dvaria n ce 

— 2ndWCEvarianoe 

— 2ndWH Evaria nee 


Fig 13. Second order means response for the 
exact, WCE and WHEP. Comparison between the 
three means at x = n I 3, <x = 1 . 


Fig 14. Second order response variances for the 
exact, WCE and WHEP. Comparison between the 
three variances at x = n!3,<y = \ . 



T=5;GRV=50;<7=1;/i=2 


variance 



— Exactvarianoe 

— 2ndWGE,variani: 

— - 2ndWH Evaria n c 


FIG 15. Second order means response for the 
exact, WCE and WHEP. Comparison between the 
three means at / = 0.5, a = 1 . 


FIG 16. Second order response variances for the 
exact, WCE and WHEP. Comparison between the 
three variances at t = 0.5, a = 1 . 


From Figs (17, 18) we found that decreasing the number of GRVs leads to a deviation of the approximate solution of WCE. 
The larger the number of GRVs, the better convergence is. 
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Fig 17. Second order means response for the exact, 
WCE and WHEP. Comparison between the three 
means at x = n / 3 and GRV = 20. 


Fig 18. Second order response variances for the 
exact, WCE and WHEP. Comparison between the 
three variances at * = k / 3 and GRV = 20. 


Finally, we show the errors between the exact solution and the approximate solutions of both WCE and WHEP techniques. 
We investigate the effects of the coefficients of SADE on the error. In Fig 19 we found that, the error between the exact 
solution and the WCE approximate solution is decreased by increasing order solution of WCE. In Fig 20 we found that, as 
the diffusion coefficient // increases, the error of second order approximate solution decreases. 



FlG 19. Error of variance response for 1st and FIG 20. Error of second order approximation variance 
2nd order approximation of WCE. Comparison of WCE. Comparison between the errors for different 
between the errors at x = n / 2. value of// at x = n / 2. 


In Fig (21-22) we found that, the error of WCE approximate solution is less than the error of the WHEP approximate 
solution. We can obtain the following conclusion: the WCE technique is better than the WHEP technique. 


T=5;GRV=50;<7=0.01;p=0,5 

Errorvar 



— 

EnufWGE 

— 

EnoiWHE 


Fig 21. Second order response variance for the 
WCE and WHE. Comparison between the errors 
of variance at x = n I 2 and jU = 0.5. 


T= 5 ;G RV= 50; <7=0.01; jj= 2 

Errorvar 



— 

ErrorWCE 

— 

ErrorWHE 


Fig 22. Second order response variance for the WCE 
and WHE. Comparison between the errors of 
variance at v = n / 2 and ju = 2. 


IX. Conclusion 

In this paper, the solution of the SADE with multiplicative white noise using WCE and WHEP techniques is introduced. The 
numerical results show that the WCE technique gives the closed form mean solution while the WHEP technique gives the 
mean solution in the form of an infinite series. Both techniques give the variance of the solution process in the form of 
infinite series. 


The numerical results not only demonstrate the accuracy of the two techniques, but also show that the WCE technique is 
more efficient than the WHEP technique. The diffusion and advection coefficient and the intensity of Gaussian white noise 
play important roles in the SADE solution. 
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The convergence of the WCE depends on increasing the number of Gaussian random variables. The convergence of the 
WHEP technique depends on decreasing the perturbation parameter <T . The larger the diffusion coefficient, the more 
convergence the approximate solution is. The effect of advection coefficient is contrary to that of diffusion coefficient, i.e., 
the smaller the advection coefficient, the more convergence the approximate solution. 


Appendix A 

To prove the formula (7.1), let we have the general form of equation (2.2) with initial condition and a periodic boundary 
condition on X not for — L <X<L but for 0 < X < 2 71 , 

u t (x,t ) = au xx (x,t) + F(x,t\ u(x, 0 ) = /(x), 

du(27T,t) (10.D 

u(0,t) = u(2x,t), —--=---, 

ax ax 

for the homogenous case of equation (10.1), the solution will be in the form u(x , t ) = ^ T n (t)Q n (x) , then we have 

n 

— = -a AT, ^ + AQ = 0, (10.2) 

dt dx 

there are three cases for Ai . For A, > 0 the solution will be Q(x ) = Acos(VXx) + B sin(x/Xx) , the second case for 
Ai — 0 the solution will be Q(x) = A + Bx and the third case is Ai < 0 and the solution will be 

Q(x ) = A COsh(x/-/Lx) + B sinh (^J—Ax) , after the application of boundary conditions for the three cases we will have 
the general Fourier series of a function with 2 7T period as 


X(x) = £ X\ cos (nx) + £ X 2 n sin (nx) = £X&(x), 


(10.3) 


n =0 

with the Fourier coefficients 

In 


n =1 


n,J 


^ LK 

X 0 = — \x(x)dx. 

In 0 

f 2^- 

X^ = — [ X(x)cos(nx)Jx, 

7T ^ 

0 

^ In 

X 2 = — \ X(x) sin(nx) dx. 

IT J 


(10.4) 


Let ii(x.t) = Yu* ( t)Q* n (x) in equation (10.1) we get: Y — u J n (t)Q J n (x)-aYu* ( t)Q ] n (x)" = V F* (t)Q J n (x), 

77 dt 77 77 


n,J 


z 


n,J 


-ul (0 + a\u j n (t) - Fj 0 1 ) (x) = 0, 


and for the initial condition, 

/w=Z c ie;w, 

n,j 

u(x. o) =Y.< (°>a c’q: (°> - c - ] a: w=o. 


n,j 


n,j 


n,j 


(10.5) 


( 10 . 6 ) 
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4 - U t (0 + a K U t (0 ~ F n (0 = 0 : K ( 0 ) - C ] n = 4 , = ^ , 

at 


(10.7) 


-aA„t 


multipling both sides by e 11 and integrates then multiply by its orthognality then in the end we will get 


i(x,o=Y J 


n -0 


l 

C\e~ anlt +\e~ an2(t - s) F l n (s)ds 


cos (nx) + ^ 


n =1 


C^e~ an2t +fe~ a,,2(t ~ s> F n 2 (s)ds 


sin (nx). 
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